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OBSERVATION OF GAMMA-RAYS FROM LMC X-4 ABOVE 10 »• eV 

R.J. Protheroe and R.W. Clay 
Department of Physics, University of Adelaide 
Adelaide, South Australia 5001 

Abstract Data from the University of Adelaide air shower array 
at Buckland Park taken over a three year period have been 
analysed to search for evidence of ultra-high energy y-ray 
emission from neutron star binary X-ray sources having known 
orbital periods. We report the detection of UHE y-rays from 
LMC X-4 above 10 2 * 4 * eV, the first extragalactic object to be 
positively detected at these energies. 

1. INTRODUCTION Gamma-rays at ultra-high energies (UHE) have recently 
been observed by Samorski and Stamm (1983) and Lloyd-Evans et al. (1983) 
from Cygnus X-3 and by Protheroe et al. (1984) from Vela X-i. The 
discovery that some neutron star binary X-ray sources are emitting 
ultra-high energy y-rays has led to considerable interest and speculation 
about the origin of the y-rays (Eichler and Vestrand, 1984; Stephens and 
Verma, 1984} Hillas, 1984; Protheroe, 1984} Porter, 1984; Chanmugam and 
Brecher, 1985) and to the notion that much of the galactic cosmic radiation 
above 10 4 * eV may originate from such objects (Hillas, 1984; Wdowczyk and 
Wolfendale, 1983a, b). Following our detection of UHE y-rays from Vela X-l, 
we examined data from the Buckland Park array (Crouch et al., 1981) to test 
for an excess of air showers from regions around the galactic plane and 
galactic centre (Clay et al. 1984), and from the directions of known 
sources of y-rays at 100 MeV energies (Protheroe and Clay, 1984). 

In addition to Vela X-l, there are 14 neutron star binary X-ray 
sources in the southern celestial hemisphere which have known orbital 
periods. We have recently searched the Buckland Park data for 1979-81 for 
evidence of periodic emission of UHE y-rays from these objects (Protheroe 
and Clay, 1985) and here we summarise the results of this search. 

2. DATA ANALYSIS For each object, we examined the air showers which had 
apparent arrival directions within angle 9«= of the source direction. This 
"resolution angle", chosen to optimise the signal to noise ratio, ranged 
from 2.1° to 2.7° depending on the source declination. To test for 
evidence of periodic emission of UHE y-rays in phase with the orbital 
motion of each binary system, the arrival time of each air shower was 
converted to an orbital phase, f (0 < f < 1), after applying a heliocentric 
correction. Cosmic ray events would be expected to have a uni form phase 
distribution and the hypothesis of uniformity was tested using the 
statistic ¥„ proposed by Protheroe (1985a) specifically for testing 
circular data for uniformity while being powerful against alternatives in 
the form of a uniform distribution plus a small narrow distribution at an 
arbitrary phase. Y„ exceeded the 17. critical value only for LMC X-4. For 
the 13 other objects, 95% upper limits were obtained for the y-ray flux and 
luminosity and are given by Protheroe and Clay (1985). For LMC X-4, 63 
events were recorded within Gc^.S 0 of the source direction and, using 

the ephemeris of Kelley et al. (1983), Y«» was 7.55 compared with the 5X 
and 17. critical values of 7.32 and 7.54 respectively. 

The LMC X-4 observation may actually be more significant than 



248 


06 2.6-10 


suggested by the 1% 
significance level obtained 
using the Y„ statistic if 
the peak in the phase 
distribution is broader 
than mt" 1 of the period. 

The phase distribution of 
events from within 2.5° of 
the direction of LMC X-4 is 
plotted in Fig. 1. There 
is a significant excess of 
events between phases 0.90 
and 0.95 where we observe 
12 events when we expect 
2.65 events. The 
probability, p, of this 
occuring by chance in a 
particular bin, obtained 
using the Poisson distribution, is 2.12x10"*. The phase distribution is 
shown binned in 20 bins in Fig. 1, but was previously binned in 10 bins and 
this binning and re-binning must be taken into account when assigning a 
significance level for the observation. For LMC X-4 alone, the probability 
must then be increased to l-(l-p) <*®**®> = 6 . 63xl0 -* corresponding to a 
3.2 standard deviation observation. Since 14 sources were examined for UHE 
*- ra y emission, the probability that the observed phase distribution for 
LMC X-4 resulted from a chance fluctuation of cosmic ray background events 
is l-U-p)®©«**» = 0.009. It is therefore likely (99.17. confidence) that 
the phase distribution is not the result of a background fluctuation. 

We attribute the excess between phases 0.90 and 0.95, i.e. 9. 4+3. 5 
events, to Jf-rays. This corresponds to a time averaged integral flux of 
<4.6±1.7)xl0"*» photons m"* s“* above 10** eV, the effective threshold 
of the array for showers incident from declination -66°. Assuming an E~ a 
photon spectrum at the source and taking account of electron-photon 
cascading in the 3 K microwave background we obtain a luminosity of -vlO** 
erg s~* per decade which is of the same order of magnitude as that emitted 
at 2-11 keV energies (Bradt and McClintock, 1983). 

3. DISCUSSION Since LMC X-4 is in the Large Magellanic Cloud, some 50 kpc 
away, we would expect JT-rays below about 10** eV to be attenuated on 
traversing the 3 K microwave background. This topic is discussed in an 
accompanying paper (Protheroe, 1985b). Recent extensions to the Buckland 
Park array (Prescott at al, t 1983) have made measurements below 10** eV 
possible for this declination. 

LMC X-4 is a massive binary system with a 1.408* period. The 
inclination of the system is •*66® (Kelley et al., 1983) and the X-ray 
source is eclipsed between phase 0.92 and 0.08. If UHE protons are 
accelerated in the region of the compact object UHE Jf-rays could be 
produced by nuclear interactions in the atmosphere of the companion star. 

In this case, we would expect tf-ray emission at orbital phases -*0.92 and 
•*0.08 when our line of sight to the neutron star just grazes the star's 
surface. The observation of UHE JT-rays at a phase of 0.90-0.95 is 
consistent with this picture although the absence of UHE emission at phases 
0.05-0.10 is puzzling (a similar situation exists for Cygnus X-3 and Vela 



Fig. 1. Phase distribution of the 63 air showers 
observed within 2.5» of LMC X-4. The ephemeris 
of Kelley et al. (1983) was used. The dashed line 
shows the expected phase distribution of the 
cosmic ray background events. (Reproduced from 
Protheroe and Clay, 1985). 
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X-l where only one burst of UHE Jf-rays is observed per orbit). 



Fig. 2. The LMC X-4 primary star and inclined orbit of the neutron star 
drawn to scale based on the analysis of Kelley et el. (1983). Possible 
location of material responsible for obscuration of the primary between 
phases 0.6 and 0.9 (Huchings et ai. f 1978) is shown cross-hatched, assuming 
it lies in the neutron star orbit, for an orbital phase of (a) 0.6 and (b) 
0.9. 

The star appears not to fill its Roche lobe and to have a rather low 
stellar wind (Huchings et al. f 1978). Mass transfer may occur through a 
trailing accretion stream which may feed the accretion disc. Evidence for 
this comes from variable obscuration of the companion star between orbital 
phases 0.6 and 0.9 (Huchings et ai«, 1978) which could correspond to matter 
trailing behind the neutron star by between 0.2 and 0.3 of an orbit (see 
Fig. 2). If this matter is of considerable extent above the orbital plane 
it could be the target material for interactions of UHE protons produced 
near the neutron star. Such a scenario could produce f-rays only at a 
phase of about 0.9 as observed. 
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Fig. 3. The- statistic derived from the phase distribution of the 63 
events arriving within 2.5° of LHC X-4 plotted as a function of 
precessional period. The precessional periods found by: (a) Ilovaisky et 
el, (1984), and (b) Lang et el, (1981) are indicated. 


The X-ray source has high and low states associated with a ^30* 
period (Lang et al. f 1981) attributed to precession of an accretion disc. 
The 30* period is also seen at optical wavelengths and an ephemeris based 
on observations taken over seven years has recently been published by 
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Ilovaisky et al. (1984). We have searched for evidence of modulation by 
the 30* period. We plot in Fig. 3 Y M against precession period for 
periods in the range 28.5“ to 32.5“. Y*» displays a prominant feature 

near 30.6“, highly suggestive of modulation by the precession period, and 
exceeds the 5% critical value for some periods within l<r of that found by 
Lang et al. (1981). However, for the period found by Ilovaisky et al. 
(1984), P = 30.40- ± 0.03“, Y M is less than the 57. critical value. 

There is therefore no firm evidence for modulation of the UHE Jf-rays by the 
30“ period although such modulation may be revealed in subsequent data. 

If the UHE Jf-rays are modulated with the precession period the light curve 
may be expected to be symmetrical about the time of X-ray minimum. A 
search for such modulation has also been carried out using the ephemeris of 
Ilovaisky et al. (1984) and was negative. 

4. CONCLUSION We have searched the data taken during 1979-81 with the 
Buckland Park air shower array for evidence of UHE y-ray emission by 
neutron star binary X-ray sources. In addition to our detection of 
Vela X-l reported elsewhere (Protheroe et al., 1984), we have detected UHE 
y-ray emission from LMC X-4 modulated with the 1.408“ period. 
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